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Abstract — Images obtained with coiierent illumination, as is 
the case of sonar, ultrasound-B, laser and Synthetic Aperture 
Radar - SAR, are affected by speckle noise which reduces 
the ability to extract information from the data. Specialized 
techniques are required to deal with such imagery, which has 
been modeled by the G" distribution and under which regions 
with different degrees of roughness and mean brightness can be 
characterized by two parameters; a third parameter, the number 
of looks, is related to the overall signal-to-noise ratio. Assessing 
distances between samples is an important step in image analysis; 
they provide grounds of the separability and, therefore, of the 
performance of classification procedures. This work derives and 
compares eight stochastic distances and assesses the performance 
of hypothesis tests that employ them and maximum likelihood 
estimation. We conclude that tests based on the triangular 
distance have the closest empirical size to the theoretical one, 
while those based on the arithmetic-geometric distances have the 
best power. Since the power of tests based on the triangular 
distance is close to optimum, we conclude that the safest choice 
is using this distance for hypothesis testing, even when compared 
with classical distances as Kullback-Leibler and Bhattacharyya. 

Index Terms — image analysis, information theory, SAR im- 
agery, speckle noise, multiplicative model, contrast measures. 

I. Introduction 

SONAR, laser, Synthetic Aperture Radar (SAR) and ul- 
trasound B-scanners are examples of sensing devices that 
employ coherent illumination for imaging purposes. In general 
terms, the operation of these systems consists of sending 
electromagnetic pulses towards a target and analyzing the 
returning echo. In particular, the intensity of the echoed signal 
plays an important role, since it depends on the physical 
properties of the target surface Therefore, an accurate 
modelling of the echo intensity, as well as its associated noise, 
is determinant to set the extent of the imaging capabilities of 
a given sensing system. 

Noise is inherent to image acquisition. An important source 
of noise when coherent illumination is used is due to the 
interference of the signal backscattering by the elements of 
the target surface. As a consequence of such interference, 
the returning signal becomes contaminated with fluctuations 
on its detected intensity. These alterations can significantly 
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degrade the perceived image quality, as much as the ability of 
extracting information from the echo data. The resulting effect 
is called speckled noise |[T|. 

Modelling the probability distribution of image regions can 
be a venue for image analysis I^J. In particular, the widely 
employed multiplicative model leads to the suggestion of the 
Q'^ distribution for data obtained from coherent illumination 
systems |^3J-|7]. 

A direct statistical approach leads to the use of estimated 
parameters for data analysis, but a single scalar measure would 
be more useful when dealing with images. Such measure 
can be refereed to as "contrast" if it provides means for 
discriminating different types of targets [4^1, f5\, fSl. Suitable 
measures of contrast not only provide useful information about 
the image scene but also take part of pre-processing steps in 
several image analysis procedures f9). 

The derivation of expressive contrast measures is important 
for image understanding. This can be easily done when dealing 
with optical information, since contrast mainly depends on 
brightness. In the speckled data case, the main image feature 
is the roughness. Therefore, contrast measure should take it in 
account. Nonparametric methods and basic exponential mod- 
elling could not include roughness into their framework |10|. 
Indeed, simple contrast measures, such as the square ratio 
of the sample mean difference to the sum of the sample 
variances |4|, [51, can offer low computational cost. But, on the 
other hand, these simple measures can neither provide insight 
about the roughness nor offer any known statistical property 
that could furnish hypothesis testing procedures. 

Recent years have seen an increasing interest in adapting 
information-theoretic tools to image processing |[8|. In partic- 
ular, the concept of stochastic divergence fTTj has found appli- 
cations in areas as diverse as image classification |12|, cluster 
analysis fT3\, and multinomial goodness-of-fit tests |14|. Co- 
herent polarimetric image processing has also benefited, since 
divergence measures can furnish methods for assessing seg- 
mentation algorithms Q. In |15|, the Bhattacharyya distance 
was proposed as a means to furnish a scalar contrast measure 
for polarimetric and interferometric SAR imagery. 

The aim of this study is to advance the analysis of contrast 
identification in single channel speckled data. To accomplish 
this goal, measures of contrast for Q'^ distributed data are 
proposed and assessed. These measures are based on infor- 
mation theoretic divergences, and we identify the one that 
best separates different types of targets. This paper extends 
the results presented in |16|, where an exploratory analysis of 
these distances is presented. 

The article unfolds as follows. Section |n] presents the main 
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properties of the model. Section |lll] derives eight contrast 
measures and discusses their relationships. Section IV presents 
the main results, namely, the performance of these measures as 
features for target identification. Conclusions and future lines 
of research are presented in Section [V] Appendix |A] provides 
details about the distances derived for the tj" model. 

II. The Q° distribution for speckled data 

Unlike many classes of noise found in optical imaging, 
speckled noise is neither Gaussian nor additive |[T]. Proposed 
in the context of optical statistics, the most successful approach 
for speckle data analysis is the multiplicative model, which 
emerges from the physics of the image formation |jT7). In 
particular, this model has proven to be accurate for assessing 
the distribution of the SAR return signal yj. 

Such model assumes that each picture element is the out- 
come of a random variable Z called return, which is the 
product of two independent random variables, X and Y . While 
the random variable X models the terrain backscatter, the 
random variable Y models the speckle noise. 

Coherent imaging is able to provide complex-valued infor- 
mation in each pixel |18|, but the amplitude or the intensity 
of such return is the most common format in applications. 
Without any loss of generality, in this work, only the intensity 
format for images is examined. 

Backscatter carries all the relevant information from the 
mapped area; it depends on target physical properties as, for 
instance, moisture and relief. A suitable distribution for the 
backscatter is the reciprocal gamma law |3|, X ^ T^^{a,^), 
whose density function is given by 



7-" 

rx" exp 



r(-a) 



-a, 7, x > 0. (1) 



This parametrization is a particular case of the generalized 
inverse Gaussian distribution. 

Speckle Y is exponentially distributed with unitary mean 
in single-look intensity images p9[ ; therefore a multi-look 
procedure over L independent observations furnishes intensity 
speckle that can be described by the gamma distribution, Y ~ 
r(L,i), with density given by 



fY{y;L) 



r{L) 



y 



-iexp(-i2;), y>0,L>l. (2) 



In this work, the number of looks L is assumed known and 
constant over the whole image. A detailed account of the until 
recently largely unexplored issue of estimating L is provided 
in ||20|. 

Considering the distributions characterized by densities ([T]) 
and (|2|, and that the related random variables are independent, 
the distribution associated to Z = XY can be derived and its 
density is given by 

^^(^'"'^''-)" 7°r(-a)r(L) " (^ + ^^) ' 

-q;,7, z > 0,i > 1. (3) 

We indicate this situation as Z ^ Q''\a,'j, L). As shown 
in ||6), this distribution can be used as an universal model for 
speckled data. Since it has the gamma law as a particular case. 



homogeneous targets can be well described |[3). This model 
can also characterize extremely heterogeneous areas which 
are left unexplained by the K distribution [18J , for instance. 
Moreover, it is as effective as the K law for modelling 
heterogeneous data. A multivariate version of this distribution 
is presented in [181, and its application to image classification 
is discussed in pi| . 

The rth moment of Z is expressed by 



L 



r{-a-r)r{L + r) 



(4) 



if —r > a and infinite otherwise. 

Several methods for estimating parameters a and 7 are 



available, including bias-reduced procedures |22 |-p4), robust 
techniques p5) , |26| and algorithms for small samples p7). In 
this study, because of its optimal asymptotic properties p28|, 
maximum likelihood (ML) estimation is employed to estimate 
a and 7. 

Based on a random sample of size n, z = (zi, Z2, ■ ■ ■ , 
the likelihood function related to the G^{ct,'j, L) distribution 
is given by 



Thus, the estimators for a and 7, namely a and 7, respec- 
tively, are the solution of the following system of non-linear 
equations: 

-f + ^Er=i(7+i^.)-' = o, 

(5) 

where is the digamma function. However, the above 

system of equations does not, in general, possess a closed form 
solution, and numerical optimization methods are considered. 
We use the BFGS procedure, which is reportedly fast and 
accurate |29|, available in many platforms as, for instance Ox 
and R. 

Figure [T] presents a SAR image obtained by the E-SAR 
sensor over surroundings of Miinchen, Germany ||30|; its 
number of looks was estimated as 3.2. The area exhibits 
three distinct types of target roughness: (i) homogeneous 
(corresponding to pasture), (ii) heterogeneous (forest), and 
(iii) extremely heterogeneous (urban areas). Samples were 
selected and submitted to statistical analysis. Table |l] shows 
the estimates in each of these samples, as well as their size; 
the last column, namely the number of parts, will be explained 
later. Figures |2|a), (b), and (c) compare relative frequencies 
of samples to their associated Q'^ fitted densities for urban, 
forest, and pasture regions, respectively. The adequacy of the 
Q° law to speckled data is noteworthy. These samples will be 
used to validate our proposal in section [TV-BI 

As presented by |j7|, different SAR image regions can 
be discriminated using the estimated parameters of the Q'~' 
model. The expressiveness of the model and the separability 
of different samples is an open issue that we explore in this 
work. 
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Fig. 1. E-SAR image and selected regions. 

TABLE I 
Parameter estimates 



Regions 


a 


7 




# pixels 


# parts 


pasture- 1 


-15.702 


39259 


2670.32 


1235 


25 


pasture-2 


-12.698 


80320 


6866.13 


1216 


24 


pasture-3 


-11.304 


162292 


15750.39 


1602 


32 


forest 


-9.339 


661183 


79288.04 


1606 


32 


urban- 1 


-0.759 


148413 


oo. 


2005 


40 


urban-2 


-0.388 


110183 


oo. 


3481 


71 


urban-3 


-1.079 


55583 


703582.30 


4657 


95 



III. Measures of Distance and Contrast under the 

5° Law 

Contrast analysis often addresses the problem of quantifying 
how distinguishable two image regions are from each other 
In a sense, the need of a distance is implied. It is possible to 
understand an image as a set of regions that can be described 
by different probability laws. 

Information theoretical tools collectively known as diver- 
gence measures offer entropy based methods to statistically 
discriminate stochastic distributions |8|. Divergence measures 
were submitted to a systematic and comprehensive treatment 
in |[3T]-|33J and, as a result, the class of (/i, -divergences 
was proposed p3) . 

Let X and Y be random variables defined over the same 
probability space, equipped with densities fx{x;9i) and 
/y(a;;02), respectively, where 6i and 62 are parameter vec- 
tors. Assuming that both densities share a common support 
/ C M the {h, (/))-divergence, between fx and /y is defined 
by 



D'1{X,Y) 



fY{x;92) 



fY{x;02)dx 



(6) 



where 0: (0, c»)— >• [0, oo) is a convex function, h:{0,oo)^ 
[0, oo) is a strictly increasing function with h{0) = 0, and 
indeterminate forms are assigned value zero. 

By a judicious choice of functions h and </>, some well- 
known divergence measures arise. Table \n\ shows the selection 
of functions h and (j> that lead with distance measures over 
which the test powers and sizes were estimated for speckled 




00.5 X 10" 1.5 X 10" 2.5 X 10° 
Graylevel 




1X10' 2 • 
Graylevel 



(a) 



(b) 




Graylevel 

(c) 

Fig. 2. Relative frequencies (o) and fitted densities (— ) for (a) urban, 
(b) forest, and (c) pasture regions. 



data modeled by Q'^ law in p6) . Specifically, the follow- 
ing measures were examined: (i) the Kullback-Leibler diver- 
gence ||34J, (ii) the relative Renyi (also known as Chernoff) di- 
vergence of order /3 p5] , (36), (iii) the Bellinger distance (37], 
(iv) the Bhattacharyya distance f38l, (v) the relative Jensen- 
Shannon divergence [39j , (vi) the relative arithmetic -geometric 
divergence (40), (vii) the triangular distance (41], and (viii) the 
harmonic-mean distance (4T\. 

Often not rigorously a metric |42J, since the triangle in- 
equaUty does not necessarily holds, divergence measures are 
mathematically suitable tools in the context of comparing the 
distribution of random variables [43 J . Additionally, some of 
the divergence measures lack the symmetry property. Although 
there are numerous methods to address the symmetry prob- 
lem 1441, a simple solution is to define a new measure d'^ 
given oy 

L>';(x,y) + D'f(y,x) 



regardless whether D^{-,-) is symmetric or not. Henceforth, 
the symmetrized versions of the divergence measures are 
termed "distances". By applying the functions of Table |ll] 
into equation (|6]l, and symmetrizing the resulting divergences, 
integral formulas for the distance measures are obtained. 
For simplicity, in the list below, we suppress the explicit 
dependence on x and the support / in the notation. 

(i) The Kullback-Leibler distance: 

d^^{X,Y) = \ j[fx-fy)\0g(^l^y 

(ii) The Renyi distance of order /?: 



<(x,r) 



log 
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TABLE II 

(h, (/i)-DIVERGENCES AND RELATED FUNCTIONS (p AND h 



{h, 0)-divergence 


Hy) 


0(x) 


Kullback-Leibler 
Renyi (order P) 
Hellinger 
Bhattacharyya 

Jensen-Shannon 
Arithmetic-geometric 
Triangular 
Harmonic-mean 


y 

-phi log - 1)2/ + 1) , < 1/ < 

y/2,0< y<2 
-log{-y + l),0<y < 1 

y 
y 

y,o<y<2 

-log(-| + l), 0<y<2 


X log(x) 
p \ ,0 < /3 < 1 

(^)i°s(^) 

X+1 

(x-lf 
x + 1 



(iii) The Hellinger distance: 

dn{X,Y) = l-j 77^- 1-cxp (^-]^d^^\x,Y) 

(iv) The Bhattacharyya distance: 

d^{X, Y)^- log v/7Wv) = - log (1 - dn{X, Y)) 



(v) The Jensen-Shannon distance: 
1 
2 



dis(^,>") 



/xlog 
/y log 



2/y 



jY + fx 

(vi) The arithmetic-geometric distance: 

d^o{X, Y) = \j{fx + Iy) log ^ + 

(vii) The triangular distance: 

d,(X,Y)= j^^""-^^^' 



fx + fY 



(viii) The harmonic-mean distance: 



dHM{X,Y) = -l0g(/ 

= -logfl- 



VxfY 
fx+fY 

djjX^Y) 
2 

Alternatively, the distances can be put under the {h, (f)- 
formalism. The distances derived from symmetric divergences 
inherit the same h and (f) functions. For the remaining dis- 
tances, specifically tailored h and functions can be found as 
shown in Table Hill 

Provided that the concerned random variables follow the 
5° law with parameter vectors 6^ = (ai,7i,Li) and 6^ — 
(q;2,72, L2), particular expressions for the discussed distances 
can be achieved. After adequate considerations, the integral 
forms of some of the distances furnish closed expressions. 
Appendix [a] details the mathematical manipulations employed 
to derive the Kullback-Leibler, Renyi of order /3, Hellinger, 
and Bhattacharyya distances between two distributed ran- 
dom variables. By contrast, no corresponding closed form 
expressions were found for the triangular, Jensen-Shannon, 
arithmetic-geometric, and harmonic-mean distances. In order 
to evaluate them, numerically quadrature routines available for 
the Ox programing language were employed 145). 



When considering the distance between same distributions, 
only their parameters are relevant. In this case, parameter 
vectors 6^ and 62, replace random variables symbols X and 
Y as the arguments of divergence and distance measures. This 
notation is in agreement with that of |33|. 

Figure [3] depicts plots for the distances d^{6i,62) between 
■^0, where Bi = (ai,7i,8) and 62 = (-12,11,8). Parameter 
ai ranges in the interval [—14,-10] and 71 was selected, 
using equation Q, so that its associated distributed random 
variable has unit mean: 



7i 



LT{-a^)V{L) 
r(-ai-l)r(L + l) 



-ai — 1. 



(8) 



The obtained curves indicate that Hellinger and Bhat- 
tacharyya distances exhibit comparable behavior. Similarly, 
Kullback-Leibler, Renyi with /3 = 0.95, and triangular dis- 
tances have closely matching plots. 



o 

X 



e 




o 

X 



CM 

7 




-Oil 



Fig. 3. Distance measures between two 5" distiibuted random variables as 
a function of ai. 



Several convergence properties of the (/i, 0) -divergences 
were established in |33J. Under the regularity conditions 
discussed in |33 p. 380], if parameter vectors 6^ and 0-2 are 
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TABLE III 

(h, (/))-DISTANCES AND THEIR FUNCTIONS (f) AND h 



[h, </))-distance 


h{y) 


0(x) 


Kullback-Leibler 
Renyi (order /9) 
Jensen-Shannon 
Arithmetic-geometric 


2//2 

^Iog({/3-l)y + l), 0<s,< ^ 
2//2 

y 


2( 
xlog 

(^) 


(x — 

-/3(3: 

log( 


L)log(x) 

"^'"^0 < ,3 < 1 
^) + (^) 



equal, then, as m, ti cx), the quantity 

m + n /i'(0)(/)"(l) 

is asymptotically chi-square distributed with A/ degrees 
of freedom, where Oi = {On, . . . ,9i]\j) and 02 — 
{021, ■ ■ ■ ,02m) ^"6 the ML estimators of 9i and 02 based 
on independent samples of sizes m and n, respectively p3| . 

Thus, when considering the definition of the distances in 
terms of the h and (p functions, and applying the results on 
the convergence in distribution of the {h, 0)-measures to xit 
p3] , the lemma asserted below is proved. 

Lemma 1: Let the regularity conditions proposed in f33' 
p.380] hold. If > A e (0, 1) and 0i = 6>2, then 



m,n— ^oo 



2mn d'^{di,S2) 

m + n h'{0)(l)"{l) m,n~¥co 



> xli, 



(9) 



where " — !•" denotes convergence in distribution. 
Based on Lemma [T] statistical hypothesis tests for the null 
hypothesis 6i = 02 can be derived. In particular, the following 
statistic is considered: 

^ m + n ^ 

where v ~ 1/ {h' {0)(j)" (1)) is a constant that depends on 
the chosen distance. Table |IV] lists the values of v for each 
examined distance. We are now in position to state the 
following result. 

Proposition 1: Let m and n assume large values and 
S'^{9i,92) — s, then the null hypothesis 6i — 02 can be 
rejected at a level rj if Pr (xm — 

In terms of image analysis, this proposition offers a method 
to statistically refute the hypothesis that two samples obtained 
in different regions can be described by the same distribution. 

TABLE IV 

DISTANCES AND CONSTANTS V 



Distance 


V 


Kullback-Leibler 


1 


Renyi (order P) 




Hellinger 


4 


Bhattacharyya 


4 


Jensen-Shannon 


4 


Arithmetic-geometric 


4 


Triangular 


1 


Harmonic-mean 


2 



IV. Results and Discussion 

In order to assess the proposed contrast measures, both 
synthetic distributed data and actual SAR images were 
submitted to the statistical analysis suggested by Proposition [T] 
Two nominal levels of significance were considered, namely 
1% and 5%. These results are presented in sections |IV-A| 
and |IV-B] respectively. 

Usually SAR images are analyzed in square arrays of size 
7 X 7, 9 X 9, and 11 x 11 pixels. In a conservative way, we 
chose to work with the smallest sample size, i.e., windows of 
size 7x7 pixels, but we present a summary of results for 
larger windows, i.e., 9x9 and 11 x 11. 

A. Analysis with Simulated Data 

Although the distribution is specified by a and 7, SAR 
literature often employs the texture a and the mean fi. Since 
equation (|4| establishes that 

_ 7 r(-Q- i)r(£ + i) _ 7 

^~ L r(-a) r{L) ^ 1 + a 

both specifications are equivalent. Thus, prescribing the pa- 
rameter values of a S {—1-5, —3, —5, —8}, /i G {1, 2, 5, 10}, 
and L S {1,2,4,8}, a total of 64 statistically different image 
types will be used in the following assessment. 

The empirical size and power of the proposed test were 
sought as a means to guide the identification of the most 
adequate distance measure. To obtain the pursued empiri- 
cal data, Monte Carlo experiments under different scenarios 
were designed. Let two tj" distributed images be specified 
by the parameter vectors (aij/ii.L) and (a2,/X2,L), for 
L G {1, 2, 4, 8}. Four scenarios were considered in such a way 
that image pairs under scrutiny satisfy: (i) ai — a2,iJ,i 7^ fJ.2, 
(ii) ai 7^ a2,p.i = /i2, (iii) < 02, Mi < A^2, or 
(iv) ai < a2,Hi > ^2- 

Situation (i) corresponds to 71 7^ 72 and ai = 012. For the 
other three situations, let k = (l + ai)/(l + Q!2). Situation (ii) 
is 7i = K72 and ai ^ 0.2- Situation (iii) is 71/72 > k and 
ai < 0L2- Finally, situation (iv) is 71/72 < k, and a\ < 012- For 
the given selection of parameter values, pairwise combinations 
of the 64 image types furnished 96 different cases for each 
scenario (i) or (ii). Situations (iii) and (iv) offered 144 cases 
each. 

Situation (i) describes a tough task: discriminating two 
targets with equal mean brightness that only differ on the 
roughness. Situation (ii) models the situation where areas with 
equal roughness have different mean brightness. Situations (iii) 
and (iv) describe pairs of targets whose roughness and mean 
brightness are both different, but with different relations. 
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Images submitted to the suggested statistical test for ho- 
mogeneity must have their distribution parameters estimated. 
However, the employed ML estimators for Q'^ distributed data 
are often difficult to be evaluated due to numerical instability 
issues |24|. This problem was previously reported in |27|, and 
estimate censoring was proposed as a procedure to circumvent 
this situation. Given a sample, we compute the ML estimators 
(d, 7) defined in equation set ^ and apply censoring as 
explained below. 

Monte Carlo simulations were performed for each sce- 
nario and only those results where a S [10a,Q;/20] were 
recorded valid. Up to 5500 replications were considered and, 
as presented in the following tables in the 'Rep' column, at 
least 1343 valid replications were obtained. All computations 
were performed using the Ox programing language [29 1 ; in 
particular, the quasi-Newton method with analytical derivatives 
was used to obtain the estimates. 

In the following, we report the null rejection rates of tests 
whose statistics S'^ are based on the discussed stochastic dis- 
tances: Kullback-Leibler (^kl), Renyi of order /3 = 0.95 (5r), 
Bellinger (Sh), Bhattacharyya (Sb), Jensen-Shannon (5js), 
arithmetic-geometric (Sag), triangular (St), and harmonic- 
mean (S'hm)- Data was simulated obeying the null hypothesis 

"Ho : (ai,7i) = ("2,72) = (a*,7*)- 

Table [yjpresents the empirical sizes (rejection rates of sam- 
ples from the same distribution) of the tests at nominal levels 
1% and 5%. The changes in the value of 7* for a specific L do 
not alter significantly the rate of type I error Although changes 
of scale do not alter the distance between distributions, the 
application of the maximum likelihood estimation could raise 
concerns. Such estimation method is known (i) to be prone to 
severe numerical instabilities and (ii) to increase the estimator 
variance when a is reduced. In spite of these facts, the test 
performance was little affected. 

For smaller values of a* (homogeneous images), the em- 
pirical sizes are reduced; this is due to the fact that the 
distribution becomes progressively insensitive to changes 
of a, i.e., the relative difference between densities is more 
pronounced for the same variation of a when this texture 
parameter is larger The triangular distance presents the opti- 
mum performance regarding test size, since its Type I error is 
closest to the nominal values. The tests yielded the empirical 
size closest to the theoretical one as follows: triangular and 
harmonic-mean in all of the 64 situations, Jensen-Shannon in 
98.44%, HelHnger in 90.63%, Bhattachai-yya in 85.94%, Renyi 
in 75%, Kullback-Leibler in 73.44%, and, finally, arithmetic- 
geometric in 56.25%. The lowest of these cases are highlighted 
in boldface type in Table [V] It is noteworthy that the two most 
commonly employed distances, namely the Kullback-Leibler 
and the Bhattacharyya distances, presented poor performances 
when used as test statistics. 

The efficiency of the measures Sj, Sb, Sr, Su, and Sjs with 
regard to S'kl is another important fact, since it is common to 
use measures based in the Kullback-Leibler classic divergence. 

Tables [Vl] |VII) \Vm\ and |K] complete the analysis of 
the tests based on stochastic distances by presenting their 
empirical power, i.e., the rejection rates when samples from 
different distributions are contrasted. 



Table VI presents the empirical power of the tests at 1% and 
5% nominal levels when ai — 0,2 and /ii 7^ ^2- This situation 
evaluates the effect of the change of mean gray level while 
keeping the roughness constant. For fixed L, the test power 
increases as the ratio ^2 111 increases. Additionally, increasing 
the number of looks enhances the power of the test. The power 
is larger for smaller values of a, i.e., in homogeneous targets, 
which is in agreement with the aforementioned sensitivity 
of the distribution to the texture parameter. In general, the 
empirical power is high. For example, it is greater than 61.89% 
for L > 4. 

In summary, these tests are able to recognize images of same 
roughness with different mean brightness. The arithmetic- 
geometric distance provides the best test for small values of 



L (highlighted in boldface type in Table VI when there are no 
matching situations). 

It is noteworthy that as L increases there is a threshold 
for which all tests exhibit the same performance. The more 
homogeneous the target, the smaller this threshold is. As 
expected, it is easier to perform sound statistical tests on 
homogeneous areas than in heterogeneous or extremely het- 
erogeneous targets. More looks are needed in the latter cases 
for attaining the same power. 



Table VII presents the empirical power of tests at 1% and 
5% nominal levels for the case of equal mean brightness 
(/zi — /i2) but different roughness (ai 7^ a-i)- The test based 
on the arithmetic-geometric distance consistently has the best 
performance regarding this criterion, with a few situations 
where other tests match it. 

As previously said, the task of discriminating targets of 
same mean brightness, but different roughness is tough. There 
are situations where the power of the best test is as low as 0.71 
(when 0L\ — —5, = —8, 71 = 8, 72 = 14 and L — X) but, 
for fixed ai 7^ a-i the power increases with the rate i\lii and 
with the number of looks. The former situation is illustrated 
in Figure |4] 




2 3 
Intensity 



(b) 



Fig. 4. A tougti problem: (a) densities of tlie ^''(-5, 8, 1) andg°{-8, 14, 1) 
distributions (solid and dashes, respectively) and (b) the associated data (upper 
and lower half, resp.) 

The worst cases, i.e., those with smallest power, are related 
to small values of a which correspond to homogeneous areas. 
The distance between Q'^ distributions becomes less sensitive 
to different roughness in such targets. 
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Rejection rates of (h, < 



TABLE V 

-DIVERGENCE TESTS UNDER i^o : («l,7l) : 



("2,72) = (a*, 7*). a* S {-1.5,-3,-5,-8} 



1% nominal level 



5% nominal level 



-1.5 0.5 



-3.0 2 



-5.0 4 



40 



-8.0 7 



70 



■JKL 

1.06 
1.27 
1.39 
1.71 
0.71 
1.20 
1.59 
1.60 
0.47 
0.51 
0.92 
1.20 
0.54 
0.62 
0.77 
1.27 
1X28" 
0.57 
0.64 
1.00 
0.33 
0.40 
0.50 
0.96 
0.33 
0.13 
0.32 
0.61 
0.08 
0.36 
0.52 
0.66 



Sh St 



0.65 
0.82 
0.86 
1.16 
0.38 
0.81 
1.13 
0.95 
0.34 
0.35 
0.61 
0.81 
0.29 
0.49 
0.53 
0.89 
0.23 
0.47 
0.36 
0.75 
0.24 
0.34 
0.27 
0.67 
0.26 
0.08 
0.29 
0.42 
0.08 
0.31 
0.47 
0.51 



0.40 
0.52 
0.42 
0.58 
0.25 
0.47 
0.55 
0.51 
0.22 
0.16 
0.31 
0.33 
0.22 
0.25 
0.31 
0.44 
0.18 
0.31 
0.27 
0.37 
0.14 
0.31 
0.16 
0.37 
0.07 
0.08 
0.06 
0.13 
0.08 
0.22 
0.17 
0.31 



0.71 
0.88 
1.00 
1.27 
0.46 
0.90 
1.23 
1.11 
0.41 
0.37 
0.70 
0.92 
0.35 
0.51 
0.61 
1.01 
0.23 
0.47 
0.48 
0.85 
0.33 
0.34 
0.39 
0.77 
0.26 
0.08 
0.32 
0.48 
0.08 
0.36 
0.50 
0.53 



0.52 
0.69 
0.57 
0.87 
0.31 
0.60 
0.82 
0.73 
0.31 
0.30 
0.41 
0.57 
0.25 
0.39 
0.49 
0.70 
0.18 
0.41 
0.36 
0.63 
0.14 
0.34 
0.23 
0.55 
0.26 
0.08 
0.26 
0.24 
0.08 
0.27 
0.32 
0.44 



OHM 

0.50 
0.60 
0.53 
0.80 
0.29 
0.56 
0.75 
0.64 
0.31 
0.23 
0.41 
0.46 
0.29 
0.37 
0.43 
0.53 
0.18 
0.41 
0.34 
0.53 
0.14 
0.34 
0.23 
0.49 
0.26 
0.08 
0.23 
0.24 
0.08 
0.27 
0.35 
0.42 



Sag Sr 



1.79 
2.34 
2.14 
2.57 
1.38 
1.99 
2.49 
2.47 
0.88 
0.83 
1.39 
1.86 
0.80 
0.97 
1.39 
2.07 
0.64 
0.66 
0.84 
1.58 
0.57 
0.83 
0.78 
1.39 
0.46 
0.29 
0.47 
0.87 
0.23 
0.54 
0.73 
1.00 



0.98 
1.10 
1.30 
1.67 
0.58 
l.ll 
1.54 
1.55 
0.47 
0.49 
0.88 
1.09 
0.48 
0.55 
0.71 
1.18 
1128" 
0.47 
0.61 
1.00 
0.33 
0.34 
0.43 
0.90 
0.33 
0.13 
0.32 
0.59 
0.08 
0.36 
0.52 
0.64 



Skl Sh 



458 
5.97 
6.03 
6.51 
4.22 
5.97 
6.35 
6.84 
3.01 
3.21 
4.67 
5.78 
2.96 
3.44 
4.96 
5.64 

2.64 
2.98 
4.90 
2.05 
2.88 
2.95 
450 
1.89 
1.75 
1.90 
3.45 
1.92 
2.23 
2.77 
3.63 



3.31 
4.69 
5.28 
5.51 
3.05 
4.92 
5.40 
5.82 
2.19 
2.36 
3.97 
4.98 
2.45 
2.87 
3.88 
4.69 
1.70 
2.14 
2.46 
3.96 
1.62 
2.41 
2.38 
3.79 
1.56 
1.46 
1.69 
2.84 
1.77 
1.96 
2.50 
3.15 



St 

2.21 

3.29 

3.69 

4.04 

1.98 

3.15 

3.99 

4.48 

1.47 

1.85 

2.72 

3.41 

1.69 

2.40 

2.88 

3.39 

1.56 

1.67 

1.82 

3.15 

1.14 

1.89 

1.69 

2.51 

1.30 

0.96 

1.32 

2.12 

1.54 

1.52 

2.01 

2.19 



3.56 
4.96 
5.46 
5.90 
3.30 
5.20 
5.63 
6.20 
2.41 
2.52 
4.17 
5.19 
2.55 
3.00 
4.12 
4.89 
1.88 
2.23 
2.55 
4.25 
1.81 
2.54 
2.51 
4.02 
1.69 
1.50 
1.72 
3.06 
1.85 
2.01 
2.68 
3.32 



2.85 
410 
4.60 
4.89 
2.63 
4.06 
4.85 
5.22 
1.91 
2.17 
3.56 
4.21 
1.97 
2.61 
3.35 
4.17 
1.56 
1.95 
2.23 
3.68 
1.38 
2.20 
2.13 
3.32 
1.50 
1.25 
1.61 
2.49 
1.69 
1.70 
2.39 
2.77 



^HM Sag Sr 



2.54 
3.65 
4.15 
4.46 
2.38 
3.57 
4.39 
4.86 
1.82 
1.94 
3.25 
3.78 
1.94 
2.52 
3.24 
3.82 
1.56 
1.92 
2.02 
3.47 
1.24 
2.10 
1.90 
2.85 
1.30 
1.25 
1.49 
2.38 
1.54 
1.70 
2.27 
2.64 



6.27 
7.44 
7.69 
7.77 
5.91 
7.98 
7.75 
8.17 

"TeT 

4.23 
5.97 
7.07 
3.76 
4.67 
6.34 
6.84 
2.75 
3.12 
4.05 
6.30 
2.67 
3.46 
423 
5.95 
2.54 
2.05 
2.72 
433 
2.39 
2.90 
3.44 
472 



4.37 
5.76 
5.81 
6.39 
4.01 
5.76 
6.16 
6.77 

"T88" 
3.14 
4.56 
5.57 
2.77 
3.30 
4.71 
5.52 

"2:oT 
2.55 
2.89 
4.71 
2.00 
2.81 
2.74 
4.38 
1.76 
1.67 
1.87 
3.32 
1.85 
2.10 
2.74 
3.55 



Rep 
4802 
5347 
5473 
5496 
4792 
5326 
5468 
5497 
3190 
4329 
5113 
5419 
3140 
4327 
5098 
5421 
2178 
3177 
4398 
5077 
2098 
3234 
4374 
5094 
1536 
2394 
3422 
4575 
1299 
2241 
3434 
4512 



Table VIII presents the empirical power of tests at 1% and 
5% nominal levels for the case ai > a2 and /^i < /i2- In 
general, the powers are large in this case but the best test 
regarding this criterion is the one based on the arithmetic- 
geometric distance. As expected, the power increases with the 
number of looks and with the parameter difference. 



Table IX presents the empirical powers of test at 1% and 5% 
nominal levels for the case ai > a2 and fii > ii2- Considering 
L fixed, it suggests that the empirical test power is nearly the 
same for a value of the ratio 72/71- Moreover, these empirical 
powers increase with the number of looks L. 

Table |X] illustrates the performance of the tests with respect 
to the sample size. It shows the rejection rates in the same 



situation which Table VI reports in detail for N = 49, i.e., 
L = 1, 72/71 = 2 and a G {-1.5,-3,-5,-8}. This table 
shows that when the sample size varies N e {49,81,121} 
the bigger the sample the more powerful all the tests are and, 
therefore, better discrimination is achieved. 

B. SAR Data Analysis 

In this section we use the data presented in Figure [T| and 
analyzed in Table [l] as a means for validating the simulation 
results obtained in section ITV-AI 

Each of the seven labelled regions, i.e, urban-1, -2, -3, forest 
and pasture- 1, -2, -3, was partitioned into disjoint 7x7 pixel 
samples; the number of samples (parts) is presented in the last 
column of Table U 

All pairs of parts (both from the same and different regions) 
were submitted to the proposed statistical tests. Pairs coming 
from the same region served to compute the Type I errors. 



TABLE X 

Rejection rates of {h, (^)-divergence tests under 
Hx : (ai,7i) 7^ (02, 72) WHERE oi = Q2 G {-1.5, -3, -5, -8}, 
L = I, 71/72 = 2 AND SAMPLE SIZE 







1% nominal level 


5 


% nominal level 






a 


N 




St 


Sb 


5 AG 


Skl 


St 


Sb 


Sag 


Rep 


-1.5 


49 


28.56 


22.15 


27.68 


31.29 


50.36 


45.59 


49.93 


52.95 


5133 




81 


53.24 


48.79 


52.51 


55.52 


7461 


72.23 


74.44 


76.00 


5391 




121 


75.04 


72.68 


74.84 


75.85 


90.15 


89.27 


89.98 


90.59 


5461 


-3 


49 


39.34 


30.70 


37.87 


42.63 


62.42 


57.08 


61.31 


64.20 


4143 




81 


68.03 


63.52 


67.26 


70.14 


85.30 


83.24 


84.93 


86.26 


4679 




121 


88.64 


86.89 


88.42 


89.42 


96.89 


96.29 


96.83 


97.15 


5010 


-5 


49 


4488 


35.28 


43.24 


49.31 


68.28 


62.53 


66.88 


70.44 


3427 




81 


75.65 


70.50 


75.01 


77.89 


91.00 


88.70 


90.49 


91.59 


3922 




121 


9430 


92.73 


93.93 


94.58 


98.39 


98.13 


98.36 


98.55 


4280 


-8 


49 


48.62 


38.61 


46.83 


53.17 


72.02 


65.40 


70.87 


74.26 


2945 




81 


80.10 


74.80 


78.90 


82.15 


92.74 


90.63 


92.31 


93.29 


3266 




121 


96.00 


94.77 


95.81 


96.55 


99.21 


99.01 


99.15 


99.45 


3651 



while pairs extracted from different regions were used to 
calculate the Type II errors under the hypothesis ai > a2 and 
A*i > /^2, in accordance with the estimates shown in Table |l] 

Table XI presents the observed rejections rates of samples 
from the same region. The results show that all the tests 
have excellent performance for pasture and forest regions with 
respect to this criterion. In urban scenarios the results show 
that ^KL, St and 5js maintain the good performance. Two 
additional observation are noteworthy: the Sag test shows 
optimal size only at the 1% level, and the ^h, ^b, and Sr 
classical test present an instability in the estimated size, as 
well as Sac when the nominal level is 5%. 

As previously mentioned, the Q'^ distribution is quite sen- 
sitive to the roughness parameter in extremely heterogeneous 
situations, and small random fluctuations may produce test 
statistics leading to rejection. The test size decreases with the 
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TABLE VI 

Rejection rates of (h, </))-divergence tests under Hi (ai,7i) ^ (02, 72), 01 = 02 



■ a 6 {-1.5,-3,-5,-8} 



1% nominal level 



5% nominal level 



-1.5 



a L 

u 

2 1 
2 
4 
8 
1 
2 
4 



2.5 



5.0 



S'kl 
28.56 
48.90 
68.25 
79.48 
53.07 
79.50 
93.20 
97.16 
98.64 
99.98 
100.00 
100.00 



& 

26.46 
47.10 
66.43 
77.61 
51.02 
78.46 
92.65 
96.87 
98.47 
99.98 
100.00 
100.00 



St 
22.15 
42.65 
61.89 
73.14 
46.66 
75.19 
90.59 
95.65 
98.04 
99.98 



Sb 
27.68 
48.66 
67.91 
78.96 
52.46 
79.28 
93.09 
97.13 
98.64 
99.98 



S,s 
24.59 
45.51 
64.90 
76.06 
49.50 
77.46 
92.03 
96.34 
98.33 
99.98 



Sum 
24M 
45.56 
6453 
75.61 
49.89 
77.56 
92.05 
96.20 
98.37 
99.98 



100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 



Sag 
31.29 
51.97 
71.11 
82.18 
56.36 
81.33 
94.20 
97.85 
98.85 
99.98 
100.00 
100.00 



Sr_ 

28.37 
48.88 
68.16 
79.37 
52.99 
79.50 
93.18 
97.14 
98.62 
99.98 
100.00 
100.00 



S'kl 
50.36 
71.93 
85.49 
92.29 
75.51 
92.35 
98.12 
99.40 
99.83 
100.00 
100.00 
100.00 



Sh^ 

48.88 
71.04 
8478 
91.71 
74.23 
91.94 
98.03 
99.31 
99.81 
100.00 
100.00 
100.00 



St 
45.59 
68.35 
83.03 
90.14 
71.79 
90.81 
97.70 
99.11 
99.71 



Sb 
49.93 
71.56 
85.31 
92.09 
75.02 
92.18 
98.10 
99.36 
99.81 



Sis 

Aim 

70.08 
84.16 
91.11 
73.33 
91.53 
97.92 
99.24 
99.77 



Shm 
W3b 
69.91 
83.98 
90.91 
73.21 
91.48 
97.89 
99.20 
99.75 



100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 



Sag 
52.95 
73.23 
86.77 
93.03 
76.95 
92.90 
98.23 
99.47 
99.84 
100.00 
100.00 
100.00 



Sr 

50.24 
71.85 
85.46 
92.27 
75.38 
92.27 
98.12 
99.36 
99.83 
100.00 
100.00 
100.00 



Rep 
5133 
5397 
5487 
5498 
5133 
5409 
5486 
5496 
5152 
5397 
5486 
5496 



-3.0 



2 1 
2 

4 
8 

2.5 1 
2 

4 



5.0 



39.34 
68.87 
89.82 
97.41 
70.15 
93.78 
99.43 
99.96 
99.81 
100.00 
100.00 
100.00 



36.33 
67.19 
89.12 
97.06 
67.60 
93.37 
99.40 
99.96 
99.78 
100.00 
100.00 
100.00 



30.70 
63.25 
87.08 
95.89 
62.23 
91.97 
99.17 
99.95 
99.71 
100.00 
100.00 
100.00 



37.87 
68.46 
89.74 
97.34 
69.00 
93.72 
99.43 
99.96 
99.78 
100.00 
100.00 
100.00 



33.91 
65.77 
88.40 
96.59 
65.63 
92.98 
99.30 
99.96 
99.76 
100.00 
100.00 
100.00 



33.77 
66.00 
88.35 
96.61 
65.15 
92.98 
99.32 
99.96 
99.76 
100.00 
100.00 
100.00 



42.63 
71.31 
91.03 
97.78 
73.33 
94.65 
99.51 
99.98 
99.88 
100.00 
100.00 
100.00 



39.13 
68.78 
89.82 
97.38 
69.83 
93.76 
99.43 
99.96 
99.81 
100.00 
100.00 
100.00 



62.42 
86.99 
97.07 
99.47 
87.74 
98.25 
99.94 
100.00 
100.00 
100.00 
100.00 
100.00 



60.54 
86.45 
96.92 
99.41 
86.94 
98.07 
99.92 
100.00 
100.00 
100.00 
100.00 
100.00 



57.08 
8451 
96.24 
99.30 
84.78 
97.67 
99.92 
100.00 
100.00 
100.00 
100.00 
100.00 



61.31 
86.88 
97.00 
99.47 
87.33 
98.23 
99.94 
100.00 
100.00 
100.00 
100.00 
100.00 



59.35 
85.61 
96.73 
99.39 
86.31 
97.96 
99.92 
100.00 
100.00 
100.00 
100.00 
100.00 



58.87 
85.51 
96.71 
99.38 
85.85 
97.90 
99.92 
100.00 
100.00 
100.00 
100.00 
100.00 



64.20 
88.01 
97.49 
99.50 
88.79 
98.35 
99.94 
100.00 
100.00 
100.00 
100.00 
100.00 



62.18 
86.90 
97.05 
99.47 
87.65 
98.25 
99.94 
100.00 
100.00 
100.00 
100.00 
100.00 



4143 
4879 
5294 
5451 
4120 
4858 
5294 
5451 
4142 
4868 
5280 
5451 



-5.0 



2.5 



5.0 1 
2 

4 
8 



44.88 
79.35 
96.95 
99.64 
76.25 
97.65 
99.98 
100.00 
99.94 
100.00 
100.00 
100.00 



41.35 
77.51 
96.56 
99.60 
74.33 
97.26 
99.98 
100.00 
99.94 
100.00 
100.00 
100.00 



35.28 
73.65 
95.39 
99.53 
68.30 
96.48 
99.96 
100.00 
99.91 
100.00 
100.00 
100.00 



43.24 
78.72 
96.82 
99.62 
75.20 
97.50 
99.98 
100.00 
99.94 
100.00 
100.00 
100.00 



39.04 
76.06 
96.09 
99.55 
72.04 
97.04 
99.96 
100.00 
99.91 
100.00 
100.00 
100.00 



38.40 
75.98 
96.11 
99.57 
71.32 
96.99 
99.96 
100.00 
99.94 
100.00 
100.00 
100.00 



49.31 
82.24 
97.32 
99.74 
79.03 
98.03 
99.98 
100.00 
99.94 
100.00 
100.00 
100.00 



44.67 
79.26 
96.97 
99.62 
76.07 
97.65 
99.98 
100.00 
99.94 
100.00 
100.00 
100.00 



68.28 
92.75 
99.28 
99.92 
91.24 
99.39 
100.00 
100.00 
99.97 
100.00 
100.00 
100.00 



66.09 
92.36 
99.22 
99.92 
90.66 
99.39 
100.00 
100.00 
99.97 
100.00 
100.00 
100.00 



62.53 
91.05 
99.14 
99.89 
88.52 
99.30 
100.00 
100.00 
99.97 
100.00 
100.00 
100.00 



66.88 
92.58 
99.26 
99.92 
90.95 
99.39 
100.00 
100.00 
99.97 
100.00 
100.00 
100.00 



64.78 
91.87 
99.20 
99.92 
89.94 
99.37 
100.00 
100.00 
99.97 
100.00 
100.00 
100.00 



63.99 
91.65 
99.20 
99.92 
89.56 
99.39 
100.00 
100.00 
99.97 
100.00 
100.00 
100.00 



70.44 
93.43 
99.30 
99.94 
92.31 
99.49 
100.00 
100.00 
99.97 
100.00 
100.00 
100.00 



68.02 
92.72 
99.28 
99.92 
91.21 
99.39 
100.00 
100.00 
99.97 
100.00 
100.00 
100.00 



3427 
4122 
4880 
5291 
3448 
4122 
4885 
5291 
3427 
4177 
4880 
5297 



-8.0 



2.5 



5.0 



48.62 
82.03 
98.96 
99.98 
82.44 
93.47 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



44.99 
80.58 
98.82 
99.98 
79.54 
93.30 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



38.61 
76.09 
98.64 
99.98 
73.06 
92.83 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



46.83 
81.37 
98.94 
99.98 
80.79 
93.36 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



42.51 
79.10 
98.73 
99.98 
77.44 
93.11 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



41.83 
78.39 
98.75 
99.98 
76.41 
93.11 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



53.17 
83.65 
99.22 
100.00 
85.31 
93.61 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



48.25 
81.97 
98.96 
99.98 
82.06 
93.44 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



72.02 
90.57 
99.77 
100.00 
94.27 
94.06 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



70.05 
90.12 
99.72 
100.00 
93.27 
94.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



65.40 
89.24 
99.65 
100.00 
91.82 
94.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



70.87 
90.35 
99.77 
100.00 
93.72 
94.03 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



68.15 
89.75 
99.68 
100.00 
92.65 
94.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



67.10 
89.64 
99.70 
100.00 
92.31 
94.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



74.26 
91.00 

99.77 
100.00 
94.93 
94.12 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



71.92 
90.52 
99.77 
100.00 
94.14 
94.06 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 



2945 
3522 
4336 
4964 
2899 
3569 
4336 
5005 
2706 
3543 
4336 
4926 



area roughness value and the statistic based on the triangular 
distance Sj assumes the lowest values. 



Table XII shows the tests power at 1% and 5% nominal 



levels. Table IX leads to the conclusion that the pasture areas 



are the easiest ones to differentiate from other types of land 
cover, since the power is usually highest when constrasting 
them. In both tables, the last column indicates the number of 
image parts, which is limited to those situations where feasible 
estimates were obtained in both samples. 

V. Conclusions 

This paper presented eight statistical tests based on stochas- 
tic distances for contrast identification through the variation of 
parameters a and 7 in speckled data modelled by the dis- 
tribution. Our methodology differs from previous approaches, 
since it relies on the symmetrization of the {h, -divergence 
obtained for the model. Following this approach, it was also 
possible to find compact formulas for the Kullback-Leibler, 
Bellinger, Bhattacharyya, and Renyi contrast measures. 



We presented evidence suggesting that the measures Si, 
5b, Sr, ^h, and 5js have empirical Type I errors smaller 
than the ones based on the Kullback-Leibler distance, S'kl, 
which deserves lots of attention due to its linking with the log- 
likelihood function |46|. Regarding the power of the associated 
tests, the 5ag measure presented the best performance. 

However, we observed that for a given number of looks, the 
test power performance of the proposed measures was roughly 
the same, suggesting the test based on the triangular contrast 
measure as the better tool for heterogeneity identification. Both 
synthetic and actual data analysis support this conclusion. 

The Q'^ distribution is quite appropriate for describing 
situations of extreme roughness, i.e., with values of a close 
to zero. In this situation, despite the variability, the tests were 
also efficient. Furthermore, the power, in general, improves 
with the increase in the number of looks; that is, the measures 
of contrast perform better in images with better signal-to-noise 
ratio. 

The results we presented should lead to an informed choice 
of distances in applications as, for instance, feature selection. 
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Rejection rates of {h, 



TABLE VII 

-DIVERGENCE TESTS UNDER Hi : {ai,-fi) ^ («2, 72), WITH fli 



1% nominal level 



5% nominal level 



a2 



-5.0 



-8.0 14 



16.15 
43.20 
77.18 
92.41 
29.40 
73.70 
97.11 
99.91 
37.83 
83.25 
99.61 
100.00 



Sh_ 

14.13 
41.03 
75.01 
91.03 
25.88 
70.68 
96.55 
99.89 
32.63 
81.24 
99.42 

100.00 



St 
11.64 
35.79 
69.38 
87.65 
21.80 
65.85 
95.12 
99.77 
28.18 
77.66 
99.14 



Sb 
15.14 
42.49 
76.43 
91.71 
27.09 
72.19 
96.86 
99.89 
33.92 
82.18 
99.49 



13.09 
39.27 
73.08 
89.69 
24.35 
69.02 
95.94 
99.87 
30.31 
79.83 
99.35 



1X50 
39.46 
72.54 
88.96 
24.95 
68.92 
95.69 
99.81 
31.09 
79.69 
99.21 



100.00 100.00 100.00 100.00 



Sag 
18.56 
46.60 
79.56 
94.06 
33.07 
77.18 
97.62 
99.96 
43.24 
85.82 
99.70 
100.00 



Sr_ 

15.89 
42.97 
77.01 
92.28 
28.56 
73.36 
97.03 
99.91 
36.57 
82.99 
99.58 

100.00 



Skl 
35.38 
67.81 
90.47 
97.76 
54.13 
89.77 
99.47 

100.00 
64.71 
92.40 
99.91 

100.00 



Sh_ 

33.18 
66.41 
89.88 
97.50 
51.42 
88.68 
99.36 

100.00 
61.18 
91.92 
99.91 

100.00 



St 
30.66 
63.39 
88.54 
96.64 
47.52 
86.97 
99.18 



Sb 

67.20 
90.17 
97.59 
52.16 
89.06 
99.41 



S,s 
32.04 
65.36 
89.56 
97.26 
49.92 
88.04 
99.30 



Shm 

3Mr 

65.11 
89.39 
97.08 
49.61 
87.85 
99.28 



100.00 100.00 100.00 100.00 
57.16 62.25 59.37 58.93 
90.68 92.20 91.36 91.24 
99.88 99.91 99.91 99.91 

100.00 100.00 100.00 100.00 



Sag 
38.08 
69.53 
91.39 
98.00 
57.74 
90.79 
99.57 

100.00 
67.66 
92.80 
99.91 

100.00 



Sr_ 

34.86 
67.75 
90.43 
97.72 
53.85 
89.67 
99.45 

100.00 
64.16 
92.40 
99.91 

100.00 



-3.0 -8.0 3.5 



-5.0 



0.58 
1.70 
4.81 
12.14 
0.90 
4.41 
16.75 
45.07 



0.50 
1.04 
3.76 
10.19 
0.68 
2.80 
12.72 
39.95 



0.39 
0.83 
2.44 
6.36 
0.41 
1.83 
9.36 
30.57 



0.54 
1.14 
4.07 
10.70 
0.68 
3.09 
13.87 
41.46 



0.42 
0.99 
3.19 
8.67 
0.50 
2.40 
11.27 
36.28 



0.46 
0.99 
3.05 
7.87 
0.50 
2.46 
10.79 
34.24 



1.00 
2.66 
6.49 
15.60 
1.85 
6.18 
21.59 
51.70 



0.58 
1.54 
4.56 
11.73 
0.77 
3.94 
16.15 
44.34 



3.24 
7.14 
15.82 
30.72 
5.40 
14.34 
38.03 
69.39 



2.47 
5.62 
14.09 
28.10 
4.01 
11.57 
34.21 
66.44 



1.93 
4.10 
11.07 
23.47 
3.29 
9.30 
28.38 
61.00 



2.70 
5.86 
14.58 
28.73 
4.28 
12.01 
34.93 
67.37 



2.27 
4.90 
12.58 
26.17 
3.60 
10.43 
31.69 
64.63 



2.27 
4.45 
11.99 
24.76 
3.69 
10.12 
30.40 
62.98 



4.51 
9.37 
18.27 
34.88 
8.01 
17.68 
42.39 
72.68 



2.93 
6.66 
15.54 
30.44 
5.13 
13.74 
37.36 
68.93 



-5.0 -8.0 1.75 



0.11 
0.67 
1.17 
3.54 



0.11 
0.41 
0.75 
2.67 



0.05 
0.26 
0.39 
1.32 



0.16 
0.49 
0.86 
2.92 



0.11 
0.41 
0.62 
1.99 



0.11 
0.41 
0.52 
1.74 



0.71 
0.90 
2.18 
5.21 



0.11 
0.64 
1. 12 
3.39 



1.74 
3.55 
5.97 
12.52 



1.25 
2.95 
4.67 
10.61 



0.98 
2.28 
3.29 
7.66 



1.30 
3.06 
5.01 
10.92 



1.09 
2.65 
4.10 
9.27 



1.03 
2.43 
3.92 
8.38 



2.55 
4.52 
7.37 
15.10 



1.68 
3.40 
5.65 
12.19 



image classification, edge detection, and target identification. 
In particular, the triangular distance St is the best choice and 
its good properties do not impose an extreme computational 
burden: two ML estimates obtained by the BFGS algorithm 
and a numerical integration of a single-valued function. Using 
the Ox programming language, version 4.1, on an Intel Pen- 
tium © IV CPU at 3.20 GHz, running Windows XP, the com- 
putational time for evaluating the triangular distance between 
given samples took typically less than one millisecond. We 
used the MWC_52 pseudorandom number generator (George 
Marsaglia multiply-with-carry with the use of 52 bits), which 
has a period of approximately 2^'^^^. 

Improved estimators (bias reduction by numerical and an- 
alytical approaches, and robust versions) for the parameters 
of the Q'^ family are available (see, for instance, |22|-|26|), 



and the impact of using such estimates on the aforementioned 
distances and tests is a future line of research along with their 
extension to polarimetric distributions. 

Appendix 

Consider two random variables distributed according to 
the tj" law with parameter vectors 6i = (ai,7i,Li) and 
^2 = {oi2, 12,1^2), respectively. In this case, Kullback-Leibler, 
Renyi of order f3, Hellinger, and Bhattacharyya distances can 
be manipulated into expressions that encompass integral terms 
that are suitable for contemporary symbolic mathematical soft- 
ware |47|. The above distances are detailed and the involved 
integrals are given in closed formulas in Figure |5] Given the 
law parameter space. Figure |6] provides integral identities 
needed to derive our results. 
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TABLE Vni 

(ai, 7l) ¥= (a2, 72), WITH //I < ^l2 AND VARYING PAIRS OF a WHERE ai > a2 



10 



cti ce2 ^ L 


i'A nominal le\el 


nominal level 


Rep 


Skl Sa St Sb Sjs Sum Sag Sr 


Skl Sh St Sb Sjs Sum Sag Sr 


-3.0 -5.0 4 1 

2 

4 
8 

10 1 
2 
4 
8 


52.51 50.31 45.57 51.89 48.84 49.31 55.28 52.39 
87.54 86.83 84.43 87.54 86.04 86.45 88.77 87.54 
98.89 98.86 98.70 98.93 98.82 98.84 99.10 98.89 
99.92 99.92 99.90 99.92 99.92 99.92 99.96 99.92 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


75.02 74.06 71.74 74.60 73.32 73.36 76.68 74.83 
96.45 96.29 95.74 96.42 96.26 96.18 96.75 96.45 
99.87 99.85 99.81 99.87 99.83 99.83 99.87 99.87 
99.98 99.98 99.98 99.98 99.98 99.98 99.98 99.98 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


2594 

3661 
4757 
5250 
2572 
3734 
4738 
5263 


-3.0 -8.0 7 1 
2 

4 
8 

17.5 1 
2 

4 
8 


58.83 57.11 53.03 58.29 56.02 56.16 61.28 58.83 
93.03 92.52 91.27 93.07 91.97 92.33 93.77 93.10 
99.93 99.88 99.78 99.93 99.86 99.86 99.95 99.93 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 99.91 99.91 100.00 99.91 99.91 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


79.89 78.85 76.72 79.71 78.17 78.13 81.48 79.89 
98.17 98.01 97.59 98.17 97.91 97.88 98.20 98.17 
99.98 99.98 99.98 99.98 99.98 99.98 100.00 99.98 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


2208 
3115 
4143 
4955 
2163 
3155 
4155 
4940 


-5.0 -8.0 3.5 1 
2 
4 

8 

8.75 1 

2 
4 
8 


47.26 44.76 39.22 46.44 42.75 42.53 50.35 46.93 
87.67 86.74 83.90 87.49 85.92 86.03 88.98 87.64 
99.30 99.09 98.83 99.22 99.01 99.01 99.43 99.30 

99.98 99.98 99.96 99.98 99.98 99.98 99.98 99.98 
99.95 99.95 99.95 99.95 99.95 99.95 99.95 99.95 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


72.68 70.94 67.25 71.81 69.69 69.31 74.04 72.35 
96.00 95.78 95.22 95.93 95.55 95.55 96.26 95.97 
99.90 99.90 99.87 99.90 99.90 99.90 99.90 99.90 

99.98 99.98 99.98 99.98 99.98 99.98 99.98 99.98 
100.00 100.00 99.95 100.00 99.95 99.95 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


1841 
2677 
3855 
4771 
1841 
2677 
3803 
4822 


-1.5 -3.0 8 1 
2 
4 
8 

20 1 

2 
4 
8 


89.51 88.81 86.89 89.48 88.24 88.57 90.47 89.51 

99.52 99.48 99.31 99.52 99.43 99.43 99.60 99.52 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 

100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


96.42 96.31 95.87 96.44 96.18 96.18 96.68 96.42 
99.90 99.90 99.90 99.90 99.90 99.90 99.90 99.90 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 

100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


3851 
4775 
5290 
5442 

3858 
4774 
5298 
5442 


-1.5 -5.0 16 1 
2 
4 
8 

40 1 

2 

4 
8 


95.52 95.11 94.08 95.55 94.74 95.05 95.96 95.52 
99.98 99.98 99.98 99.98 99.98 99.98 99.98 99.98 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


98.81 98.78 98.65 98.84 98.72 98.75 98.84 98.81 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


3192 
4106 
4875 
5299 
3162 
4106 
4867 
5310 


-1.5 -8.0 28 1 
2 

4 
8 

70 1 
2 
4 
8 


97.13 96.57 95.80 97.02 96.28 96.50 97.50 97.09 
100.00 99.97 99.97 100.00 99.97 99.97 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 


99.52 99.52 99.48 99.52 99.48 99.48 99.52 99.52 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00 
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Rejection rates of (h, 0)-divergence tests under Hi 



TABLE IX 

(ai,7l) 7^ (a2,72), Ml > M2 AND varying pairs OFaWHEREai > a2 



ai 



02 



0.4 



1'/ 



nominal level 



Skl 
32.15 
62.31 
83.94 
94.27 
99.81 
100.00 
100.00 
100.00 



Sh_ 

27.33 
58.02 
82.07 
93.54 
99.77 
100.00 
100.00 
100.00 



Sj^ 

20.93 
50.10 
76.35 
91.28 
99.61 
100.00 
100.00 
100.00 



Sb_ 

28.45 
59.74 
82.89 
93.94 
99.77 
100.00 
100.00 
100.00 



24.29 
55.01 
80.03 
92.78 
99.73 



Shm 

53.54 
79.34 
92.48 
99.69 



100.00 100.00 
100.00 100.00 
100.00 100.00 



Sag 
37.20 
66.95 
86.46 
95.16 
99.81 
100.00 
100.00 
100.00 



Sr_ 

31.19 
61.90 
83.75 
94.27 
99.81 
100.00 
100.00 
100.00 



nominal level 



*5kl 
56.09 
81.73 
94.14 
98.38 
100.00 
100.00 
100.00 
100.00 



Sh^ 

52.70 
79.81 
93.67 
98.13 
100.00 
100.00 
100.00 
100.00 



St^ 

47.38 
75.74 
92.18 
97.62 
100.00 
100.00 
100.00 
100.00 



Sis 
Wf} 
78.07 
92.96 
97.98 
100.00 100.00 100.00 
100.00 100.00 100.00 
100.00 100.00 100.00 
100.00 100.00 100.00 



Sb 
53.39 

80.39 
93.99 
98.27 



'S'hm 
48.96 
77.16 

92.79 
97.87 



Sag 
61.03 
84.13 
95.14 
98.48 
100.00 
100.00 
100.00 
100.00 



Sr_ 

55.40 
81.43 
94.11 
98.38 
100.00 
100.00 
100.00 
100.00 



Re 

"2594^ 
3661 
4757 
5250 
2572 
3734 
4738 
5263 



-3.0 -8.0 1.75 



0.7 



35.96 
61.99 

87.35 
96.06 
99.95 
98.57 
100.00 
100.00 



28.85 
55.99 

84.41 
94.93 
99.86 
98.57 
100.00 
100.00 



21.42 
44.43 

76.73 
92.19 
99.72 
98.57 
100.00 
100.00 



30.34 
57.37 

85.52 
95.30 
99.86 
98.57 
100.00 
100.00 



25.72 
50.79 

81.73 
94.11 
99.77 
98.57 
100.00 
100.00 



23.82 
48.03 

79.72 
93.32 
99.77 
98.57 
100.00 
100.00 



44.11 
68.83 
89.81 
97.01 

99.95 
98.57 
100.00 
100.00 



34.74 
61.03 

87.16 
95.88 
99.95 
98.57 
100.00 
100.00 



60.64 

82.66 

95.63 
98.83 
100.00 
98.57 
100.00 
100.00 



54.30 
79.74 
94.96 
98.65 
100.00 
98.57 
100.00 
100.00 



46.47 
73.10 
93.15 
98.16 
100.00 
98.57 
100.00 
100.00 



55.25 
80.55 

95.12 
98.69 
100.00 
98.57 
100.00 
100.00 



51.27 
77.75 

94.28 
98.49 
100.00 
98.57 
100.00 
100.00 



48.37 
74.86 

93.60 
98.35 
100.00 
98.57 
100.00 
100.00 



65.94 
84.82 
96.45 
99.13 

100.00 
98.57 
100.00 
100.00 



59.42 
82.28 

95.58 
98.81 
100.00 
98.57 
100.00 
100.00 



2208 
3115 

4143 
4955 
2163 
3155 
4155 
4940 



-5.0 



-8.0 0.875 1 
2 
4 



0.35 



-1.5 -3.0 



43.02 
79.75 
97.43 

99.85 
100.00 
100.00 
100.00 
100.00 
7.12 
13.68 
22.19 
29.18 
89.53 
98.62 
99.79 
99.93 



39.16 
77.29 
96.89 

99.81 
100.00 
100.00 
100.00 
100.00 
3.17 
7.77 
15.88 
22.14 
87.27 
98.22 
99.73 
99.93 



32.16 
71.31 
95.38 
99.71 
99.95 
100.00 
100.00 
100.00 
1.17 
2.49 
7.69 
12.68 
82.24 
97.07 
99.66 
99.93 



40.41 
78.52 
97.25 

99.85 
100.00 
100.00 
100.00 
100.00 
3.53 
8.21 
16.65 
23.52 
87.95 
98.43 
99.79 
99.93 



36.34 
74.75 
96.34 

99.75 
100.00 
100.00 
100.00 
100.00 
1.87 
4.77 
12.00 
18.21 
85.23 
97.65 
99.70 
99.93 



35.47 
74.11 
96.08 

99.71 
100.00 
100.00 
100.00 
100.00 
1.51 
3.25 
9.36 
15.03 
84.11 
97.49 
99.70 
99.93 



48.13 
82.78 
97.90 

99.87 
100.00 
100.00 
100.00 
100.00 
11.94 
20.13 
28.77 
36.20 
91.89 
98.93 
99.83 
99.96 



42.53 
79.49 
97.38 

99.85 
100.00 
100.00 
100.00 
100.00 
6.21 
12.50 
21.17 
28.10 
89.14 
98.53 
99.81 
99.93 



68.17 
91.93 
99.48 

99.96 
100.00 
100.00 
100.00 
100.00 
18.54 
29.47 
41.72 
50.46 
96.86 
99.64 
99.98 
100.00 



64.86 
91.00 
99.43 

99.96 
100.00 
100.00 
100.00 
100.00 
12.96 
23.96 
36.56 
46.88 
96.32 
99.64 
99.98 
100.00 



60.62 
88.98 
99.30 

99.96 
100.00 
100.00 
100.00 
100.00 
7.87 
15.37 
28.13 
38.39 
94.89 
99.43 
99.96 
100.00 



65.83 
91.41 
99.43 

99.96 
100.00 
100.00 
100.00 
100.00 
13.66 
24.77 
37.49 
47.52 
96.63 
99.64 
99.98 
100.00 



62.79 
90.40 
99.40 
99.96 
100.00 
100.00 
100.00 
100.00 
10.41 
20.15 
33.04 
43.62 
95.88 
99.58 
99.96 
100.00 



61.71 
89.80 
99.35 

99.96 
100.00 
100.00 
100.00 
100.00 
8.49 
16.75 
29.79 
40.30 
95.31 
99.50 
99.96 
100.00 



71.43 
92.87 
99.56 

99.96 
100.00 
100.00 
100.00 
100.00 
24.51 
34.66 
47.39 
55.03 
97.51 
99.77 
99.98 
100.00 



67.90 
91.86 
99.46 

99.96 
100.00 
100.00 
100.00 
100.00 
17.74 
28.59 
40.91 
50.02 
96.81 
99.64 
99.98 
100.00 



1841 
2677 
3855 
4771 
1841 
2677 
3803 
4822 
TSST 
4775 
5290 
5442 
3858 
4775 
5259 
5442 



-1.5 -5.0 



11.15 
27.47 
52.14 
74.90 
90.83 
98.76 
99.92 
100.00 



3.67 
14.52 
41.74 
68.22 
86.75 
98.15 
99.88 
100.00 



1.20 
4.53 
25.44 
55.05 
78.27 
96.37 
99.55 
99.98 



4.07 
15.64 
43.30 
69.60 
87.63 
98.34 
99.88 
100.00 



2.13 
9.23 
35.06 
63.94 
83.74 
97.59 
99.75 
100.00 



1.25 
5.97 
28.92 
58.97 
80.61 
97.03 
99.67 
100.00 



19.42 
37.87 
60.92 
80.11 
93.55 
99.07 
99.94 
100.00 



8.90 
25.26 
50.71 
73.94 
90.32 
98.68 
99.92 
100.00 



26.66 
48.78 
71.84 
88.58 
97.50 
99.78 
99.98 
100.00 



16.32 
38.89 
67.26 
86.43 
96.74 
99.68 
99.98 
100.00 



8.00 
26.25 
58.58 
81.64 
94.24 
99.49 
99.98 
100.00 



16.92 
40.09 
67.84 
86.96 
96.90 
99.78 
99.98 
100.00 



11.75 
33.44 
63.77 
84.71 
96.05 
99.63 
99.98 
100.00 



8.27 
27.91 
60.29 
83.09 
94.94 
99.56 
99.98 
100.00 



34.12 
55.99 
76.59 
90.96 
98.17 
99.85 
99.98 
100.00 



24.81 
47.27 
71.08 
88.28 
97.41 
99.81 
99.98 
100.00 



3192 
4106 
4875 
5299 
3162 
4106 
4867 
5310 



-1.5 -8.0 



7 1 
2 

4 
8 

2.8 1 
2 
4 



17.42 
41.54 
76.84 
95.61 
93.08 
95.85 
99.93 
100.00 



4.83 
25.40 
67.42 
93.43 
88.62 
95.37 
99.91 
100.00 



1.36 
9.85 
50.92 
87.92 
78.31 
93.83 
99.84 
100.00 



5.30 
26.76 
68.87 
93.87 
89.76 
95.42 
99.93 
100.00 



2.84 
17.11 
60.20 
91.40 
84.57 
94.91 
99.84 
100.00 



1.92 
12.43 
54.84 
89.70 
80.85 
94.40 
99.84 
100.00 



27.44 
52.19 
83.87 
96.99 
95.14 
96.02 
99.95 
100.00 



14.40 
38.22 
75.74 
95.31 
92.23 
95.79 
99.93 
100.00 



3455 
61.98 
90.17 
98.68 
98.16 
96.25 
100.00 
100.00 



21.99 
52.61 
87.35 
98.34 
97.35 
96.19 
99.98 
100.00 



10.13 
39.02 
81.31 
97.66 
94.59 
96.11 
99.98 
100.00 



22.50 
53.26 
87.77 
98.46 
97.42 
96.22 
99.98 
100.00 



15.65 
46.51 
85.43 
98.14 
96.65 
96.16 
99.98 
100.00 



11.53 
41.03 
82.73 
97.86 
95.14 
96.16 
99.98 
100.00 



44.35 
68.81 
92.26 
99.16 
98.75 
96.25 
100.00 
100.00 



31.86 
60.30 
89.72 
98.64 
98.05 
96.25 
100.00 
100.00 



2715 
3524 
4291 
4990 
2715 
3518 
4308 
4995 



TABLE XI 

Rejection RATES of (/i,(/i)-divergence tests under (ai,7i) = (02,72) 





1% nominal level 


5% nominal level 1 


Regions 


'S'kl 


Sh 


St 


Sb 


Sjs 


Smh 


Sag 


Sr 


Skl 


Sh 


St 


Sb 


Sjs 


Smh 


Sag 


Sr 


P 


pasture- 1 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


300 


pasture-2 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


276 


pasture-3 


0.00 


0.00 


0.00 


0.36 


0.00 


0.00 


0.36 


0.36 


1.81 


1.81 


1.81 


1.81 


1.81 


1.81 


2.17 


2.17 


120 


forest 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.33 


0.33 


0.33 


0.33 


0.33 


0.33 


0.33 


0.33 


300 


urban- 1 


0.00 


0.13 


0.00 


0.13 


0.00 


0.00 


0.00 


75.26 


1.28 


11.15 


0.26 


11.54 


0.90 


1.54 


87.69 


87.69 


780 


urban-2 


0.00 


99.96 


0.00 


99.96 


0.00 


99.20 


0.00 


100.00 


0.48 


99.96 


0.20 


99.96 


0.28 


99.92 


100.00 


100.00 


2485 


urban-3 


4.23 


3.23 


1.97 


3.56 


2.73 


2.15 


5.64 


15.03 


12.45 


11.65 


8.96 


12.05 


10.37 


9.47 


24.75 


24.75 


4465 
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TABLE XII 

Rejection RATES OF (/l, (/))-DIVERGENCE TESTS UNDER _ffi : (ai,7i) ^ (02,72), WITH 



: M2 



Regions 
pastui'e-2 
pasture-3 
pasture-3 
pasture- 1 
;pasture-2 
: pasture-3 
: pasture- 1 
:pasture-2 
: pasture-3 
pasture- 1 
:pasture-2 
pasture-3 
pasture- 1 
;pasture-2 
: pasture-3 
■Ix forest 
■2 X forest 
■3 X forest 
X urban-2 
X urban-3 
X urban-3 



■Vis .S„ 



.''11 



pasture- 1 X 
pasture- 1 
pasture-2 
forest X 
forest X 
forest X 
urban-1 x 
urban-1 X 
urban-1 X 
urban-2 x 
urban-2 x 
urban-2 
urban-3 X 
iu-ban-3 X 
urban-3 X 
tirban- 
urban- 
urban 
urban-1 
urban-2 
urban- 1 



100.00 
100.00 
94.64 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
98.70 
98.03 
94.74 
29.68 
91.05 
98.84 



100.00 
100.00 
89.13 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
99.70 
100.00 
91.12 
93.98 
93.97 
lOO.OO 



100.00 
100.00 
87.50 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
98.10 
99.49 
80.17 
26.37 
89.92 
98.50 



100.00 
100.00 
90.48 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
99.70 
100.00 
91.79 
94.19 
94.26 
100.00 



100.00 
100.00 
88.41 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
98.50 
99.27 
87.37 
27.75 
90.34 
98.74 



100.00 
100.00 
88.41 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
99.20 
100.00 
86.06 
73.73 
92.32 
100.00 



100.00 
100.00 
96.38 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
100.00 
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(i) The KuUback-Leibler distance: 
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where fc, 



(ii) The Renyi distance of order (3 : 
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Li -1, hi = ^, Ci = \ai\ + Li for i = 1,2 and d = log {k^/ki). 
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i = 1,2. 

(iv) The Bhattacharyya distance: dB(0i, ^2) = - log {Vhh /q^ x^''^+''^y'^{hi + LixY^ (62 + L2xy^dx) . 



Fig. 5. Explicit distances under the 6" model. 
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Fig. 6. Integral identities under the 6° model. 
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